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ABSTRACT 

We present constraints on cosmological and star formation parameters 
based on combining observations of the Wilkinson Microwave Anisotropy Probe 
(WMAP) and high-redshift quasars from the Sloan Digital Sky Survey (SDSS). 
We use a semi-analytic model for reionization (Chiu and Ostriker 2000) that 
takes into account a number of important physical processes both within col- 
lapsing halos (e.g., H2 cooling) and in the intergalactic medium (e.g., H2 cooling, 
Compton cooling, and photoionization heating). We find that the Gunn- Peterson 
absorption data provide tight constraints on the power spectrum at small scales 
in a manner analogous to that derived from the cluster mass function. Assuming 
that the efficiency of producing UV photons per baryon is constant, the con- 
straint takes on the form ag^^'^ ~ 0.33 in a fiat, A-dominated universe with 
h = 0.72, n = 0.99, and Qf,h'^ = 0.024. However, the calculated optical depth 
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to electron scattering of Tes ~ 0.06 is well below the value found by WMAP of 
0.17 ± (0.04 ~ 0.07). Since the WMAP constraints on r^s are somewhat degener- 
ate with the value of the spectral index n (Spergel et al. 2003), we then permit 
the primordial spectral index n to float and fixing the best fit WMAP determi- 
nation of Qoh"^ = 0.14, while normalizing the power spectrum using WMAP. In 
addition, we allow the UV-efficiency to have time-dependence. Combining the 
WMAP constraints with the quasar transmission data, our analysis then favors 
a model with Tes = 0.111°;°^, n = 0.961°;°^, implying as = 0.831°;°^ (all at 95% 
confidence), and an effective UV-efficiency that was at least ~ lOx greater at 
z ^ 6. The implied UV-efficiencies is not unreasonable for stars, spanning the 
range from 10~^'^ ~ 10^^. These results indicate that the quasar and WMAP 
observations are consistent. If future observations confirm an optical depth to 
electron scattering r^s ~ 0.1, then it would appear that no more "exotic" sources 
of UV-photons, such as mini-quasars or AGNs, are necessary; but our analysis 
indicates that a determination of Tes > 0.17 would require a more radical solution. 

Subject headings: reionization — galaxies: formation — galaxies: quasars — 
cosmology: theory — intergalactic medium — cosmic microwave background 

1. Introduction 

The primieval spectrum of cosmic density fiuctuations on large scales is determinable 
with great precision from cosmic background radiation measurements (Page et al. 2003). In 
order to determine the total relevant spectrum, information on small scales is also necessary; 
these become nonlinear at early times, so information concerning reionization in the epochs 
20 < < 6 provide vital clues. 

Recent observations of high-redshift quasars have provided the first observational sig- 
natures of the epoch of reionization. Spectra of quasars at redshift z < 6 indicate that the 
universe was almost fully ionized up to 2; ~ 6, since even a small neutral fraction in the 
intergalactic medium (IGM) would have led to complete absorption of a quasar's continuum 
radiation. However, the first absorption spectra of quasars at higher redshift indicate that 
the abundance of neutral hydrogen increases significantly for ^ > 6. 

Reionization has generally been assumed to be caused by ionizing photons created in 
early generations of stars and/or quasars. Given this premise, hydrodynamic simulations as 
well as semi-analytic calculations seem to indicate that the process of reionization should 
occur several distinct stages. First, cosmological gas falls into deep enough potential wells 
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(caused by dark matter halos), so that they can cool and collapse to high enough densities 
to produce stars and/or quasars. The UV photons produced in this process ionizes the local 
surroundings, first within the halo itself and then outside the halo creating cosmological 
H 11 regions. This is generally referred to as the "pre-overlap" stage. As the abundance of 
these H 11 regions increases (due to additional "galaxy" formation), they eventually start to 
"overlap" so that gas in the IGM becomes exposed to multiple sources of ionizing radiation. 
After this "overlap" stage, the IGM becomes optically thin except for inside self-shielded, 
high-density clouds (those without ionizing sources). 

In this picture, the rcionization history of the universe depends on both the growth of 
density perturbations and the efficiency of star/quasar formation. The former is a complex 
function of the standard cosmological parameters — the density of the universe D,o, the 
baryon abundance fib, the expansion rate Hq, and the mass power spectral index n and 
normalization erg. The latter can be calculated with some precision based on atomic physics 
(e.g., cooling rates), but ultimately depends on some unknown parameters relating to the 
efficiency of turning mass into UV photons that can escape into the IGM. 

To this overall picture has been added the recent observations by the Wilkinson Mi- 
crowave Anisotropy Probe (WMAP) satellite (Bennett et al. 2003). In particular, the high 
values of the electron optical depth to last scattering Tcs — 0.17 ± 0.04 (Kogut et al. 2003) 
and Tcs = 0.17 ±0.07 (Spergel et al. 2003) seem to indicate a much earlier epoch or rcioniza- 
tion of Zj-ci = 17 ± 5 (Spcrgcl ct al. 2003). How can these measurements be reconciled with 
the quasar Gunn- Peterson observations? Do they indicate a source of ionizing photons that 
cannot be accounted for through standard star formation? 

In this paper, we first examine how quasar transmission measurements constrain cosmo- 
logical parameters. In particular, we use the fact that the first generation of UV-generating 
objects are in fact the tail of the distribution — the rare events that have collapsed to 
high enough density to produce stars and/or quasars. They thus provide a unique probe 
of the small scale power spectrum (at < 1 Mpc scales) in an analogous way that the X-ray 
cluster mass function probes ~ 10 Mpc scales. Our constraints are derived from quasars 
absorption measurements, using theoretical predictions based on the detailed semi-analytic 
model developed by Chin and Ostriker (2000), with some improvements derived from recent 
hydrodynamic and semi-analytical work on rcionization. 

Second, we address the question of the consistency of the WMAP and quasar mea- 
surements, paying particular attention to the degeneracy in the WMAP data between the 
optical depth t^s and the spectral index n. Combining the WMAP constraints with the 
quasar transmission data, our analysis favors a model with somewhat lower Tcs and n than 
implied by the WMAP data alone, but that are still consistent at the l-cr level. 
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The organization of this paper is as follows. In § 2, we describe the modehng of reion- 
ization and Gunn- Peterson absorption. In § 3, we compare observations from both quasars 
and WMAP to model predictions in order to constrain cosmological parameters. In § 4, we 
discuss our results and discuss their implications. We summarize and conclude in § 5. 



The details of the semianalytic model are described in Chiu and Ostriker (2000); the 
basic principles are summarized here. It is based on a two-phase model of the universe in 
which a statistical filling factor for ionized gas is self-consistently calculated. It is assumed 
that the cold, neutral phase has no sources, and evolves passively with the expansion of 
the universe. The hot, ionized phase contains the ionizing sources and evolves in line with 
local particle and energy conservation averaged over the phase. The temperature of each 
phase is calculated using standard physics, with photoheating as the source of heat in the 
ionized phase and cooling via multiple mechanisms, including H2, atomic lines, and Compton 
scattering. The time evolution is determined by particle and energy conservation. In each 
case, we only consider the regions outside of collapsed gas halos, which, of course will rise 
to the virial temperature. Such collapsed halos are calculated separately, and considered 
potential sources of ionizing radiation. 

The abundance and properties of these potential ionizing sources arc calculated on 
the basis of the Press- Schechtcr formalism (Press and Schcchtcr 1974), constrained by the 
Jeans' criterion (which utilizes the calculated gas temperatures) and by a cooling criterion 
(the cooling time must be less than the dynamical time). Cooling in the halos includes the 
important contributions from H2 cooling. Halos that statisfy the Jeans criterion and that 
can cool efficiently are sources of ionizing radiation. We calculate the luminosity of each 
ionizing source using the Schmidt law 



where /* is the fraction of baryons already turned into stars (a small correction), is the 
baryonic mass of a halo, eesc is the escape fraction from the halo, is a resolution factor 
(related to the fraction of gas that can form stars) determined through calibration (see Chiu 
and Ostriker (2000) equation (28), and below), euv is the mass to UV efficiency (where we 
have absorbed the notation in Chiu and Ostriker (2000) of Chmeuv into a single efficiency), 
and tdyn is the dynamical time of the halo. Note that this halo luminosity is somewhat 



2. Modeling Reionization and Gunn-Peterson Absorption 



2.1. Summary of the Model 
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simpler than that used in Chiu and Ostriker (2000), but is consistent with that used in 
hydrodynamic simulations. 

Calibration was done by comparing the redshift of overlap, the ionizing intensity, and the 
fraction of baryonic mass in stars between the semi-analytic results and the hydrodynamic 
simulation of Gnedin (2000ab). This simulation was chosen because it is the only published 
simulation which simulates a statistically "average" universe (as opposed to an individual 
halo collapse) with sufficient resolution to follow the overlap process, and continues at least 
redshift 4. The cosmological parameters were Qq = 0.3, Q\ = 0.7, h = 0.7, Qb = 0.04, 
n = 1, and ag = 0.9. The data used for calibration from Gnedin (2000ab) are the redshift 
of overlap ~ 7, the ionizing intensity of J21 ~ 0.3 at z = A, and the fraction of baryonic 
mass in stars /* ~ 0.04 at 2; = 4. The values of e* and euv in the semianalytic model were 
adjusted to best match these three values. The results were e* = 0.03 and euv = 1-2 x 10~^, 
which gave = 7, J21 = 0.6, and /* = 0.05. Note that Gnedin (2000ab) used e* = 0.05 and 
euv = 4 X 10~^, which are remarkably similar to the vahies in the semianalytic model given 
the differences in approach. In fact, because we take each halo as a whole, it should not be 
surprising that e* is smaller for the semianalytic model than for the hydrodynamic model. 

The determination of the escape fraction eesc bears some additional discussion. We use 
a simple Stromgren sphere approximation (see Appendix) and derive 



where fa 178 is the virial overdensity, pb is the mean baryonic density, is the recombi- 
nation rate, mn is the hydrogen mass, and Eq is 13.6 eV. In this approximation, we assume 
the halo baryons have a density run, the UV output per baryon is constant in the halo, 
the halo is in ionization equilibrium, and recombinations determine the amount of photons 
absorbed locally. 

The ionizing source luminosity defined above, then, is the luminosity seen by the IGM 
— i.e., outside of the halos. A further assumption is made in the pre-overlap stage: that at 
any particular time, the ionized volume around an isolated source is linearly related to its 
UV luminosity. The constant of proportionality is determined through the global ionization 
and energy balance, but the size of the H II regions as a function of luminosity is a linear 
relation. This assumption is consistent with previous work on cosmological H II regions (e.g., 
Shapiro and Giroux 1987). 




(2) 



where the quantity 77 is given by 




(3) 
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The outputs of the model that we use below are the fiUing factor Q, the temperature 
T4, and the ionizing intensity J21. 



2.2. Density Distribution of Cosmic Gas 

One of the most important determinants of how reionization evolves is the degree of gas 
clumping. The clumping factor, defined by 

c = (p^)/w^ (4) 

is important not only for determining the ionization balance (and hence the filling factors), 
but also for determing the thermal balance. Note that this clumping factor is only used in 
the ionized region. 

In order to determine the clumping factor, the probability density function (PDF) for 
the gas density must be known. Miralda-Escude et al. (2000) found that a good fit to the 
volume-weighted PDF as seen in hydrodynamical simulations is 



Pv{A)dA = Aexp 



(A-^/3 - Co) 
2(5o/3)2 



A-f^dA (5) 



where the over density A — Pb/Pb, = 7.61(1 + z) is related to the linear rms gas density 
fluctuation, and f3 is related to the density run at high densities. The parameters A and Co 
are found by requiring the mass and volume to be normalized to unity. 

This formula, with the parameterizations of f3 and So found in Miralda-Escudc et al. 
(2000), has been used in many others in calculating characteristics of reionization (e.g., 
Songaila and Cowie 2002, Fan et al.2002,). However, it is not often noted that the param- 
eterization is based on a single simulation with a particular set of cosmological parameters. 
In particular, Miralda-Escude et al. (2000) used the simulation reported in Miralda-Escude 
et al. (1996), which was a A-CDM model with = 0.4, fl^ = 0.6, flbh'^ = 0.015, h = 0.65, 
n — 1, and as — 0.79. Since in this paper we are considering variants in the background 
cosmology, it is certainly not sufficient to use the Miralda-Escude et al. (2000) parameter- 
ization for the gas PDF. In particular, the amount of power at small scales will depend 
significantly on all the cosmological parameters (especially erg). We therefore analyzed a sec- 
ond simulation (which uses different cosmological parameters) and to generate a second set 
of parameter fits. The second simulation is one by Cen (2002b) in which Qo = 0.3, Q\ = 0.7, 
flhh^ = 0.017, h = 0.67, n = 1, and ag = 0.9. The values of (3 and 60 are given for both 
Miralda-Escude et al. (2000) and Cen (2002b) are provided in Table 1. 
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As is clear from Table 1, the value of 5o and to a lesser extent /3 depend on cosmology. 
This is not surprising given that Sq must depend on the power spectrum. The cosmological 
dependence on P is less clear, especially since the best fit value of P depends significantly on 
the density run in collapsed halos. Because the Cen (2002b) simulation is at much higher 
resolution and includes more realistic physics, we simply fit f3 to the Cen results as a function 
of redshift, with a maximum of /3max = 2.5, corresponding to an isothermal sphere: 

/3RiMin[2.5, 3.2 -4.73/(1 + 2>4 (6) 

where wc are only considering redshifts z > 4. 

As for a prescription for finding Sq, we note that the gas PDF also predicts the fraction 
of mass in collapsed virialized halos: 

POO 

/gas PDF (collapsed) « / APv{A)dA. (7) 

The integration point Gtt^ is derived from the fact that for a singular isothermal sphere, the 
local overdensity density at the virial radius is Gtt^. Therefore, we are approximately taking 
into account all gas within virialized halos. Similarly, from linear theory, we can use the 
Press-Schechter formalism (Press and Schechter 1974) to find the same quantity if we know 
the correct filtering radius Rf: 



/ps (collapsed) J / exp 



oo 



dS, (8) 



where Sc ~ 1.69 and aRj, is the hnear RMS mass fiuctuation filtered with a tophat of radius 
Rf. However, the mass fraction should be equal by the two calculations 

/gas PDF (collapsed) = /ps (collapsed). (9) 

By analyzing both simulations, we find the following relation leads to satisfactory results: 

Rf ^ Rj/4, (10) 

where Rj is the Jeans length defined by 



^'-^ij^i^ ^^^^ 



and T is the gas temperature, p is the mean total density, is the baryons per particle. With 
this prescription for finding Sq and /3, we need only use the normalization constraints to find 
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A and Co for an arbitrary background cosmology. While by no means perfect, this procedure 
should account to first order the dependence of the gas PDF on background cosmology, and 
is certainly more accurate than a treating the evolution of Sq as independent of cosmology. 

Given the gas PDF, one can determine the clumping factor by simple integration over 
all densities. However, in reahty, there is a high density cut-off due to the fact that the very 
high density regions arc not participating in the ionization balance due to self-shielding. 
However, the detcrminination of the cutoff in our prescription is not difficult, since we have 
explicitly made a separation between the "in-halo" and "out-of-halo" calculations. Thus, we 
calculate the clumping factor by integrating 

As above, we assume that gas with overdensity > Gtt^ is within collapsed halos. 

We should note that this treatment produces smaller clumping factors (~ 5) than typ- 
ically used in semi-analytic treatments. This is because we are essentially using a "hybrid" 
between the reionization treatment of Miralda-Escude et al. (2000) and the usual clumping 
factor methods. Conceptually, Miralda-Escudc ct al. (2000) considered a density-dependent 
ionization fraction — they assumed that the universe was ionized up to a critical density, 
above which was neutral and self-shielded. Here, we are assuming that the collapsed halos 
are self-shielded, and then treat the rest of the universe using the clumping factor approach. 
As described in the appendix, at redshifts before full reionization, halos are optically thick if 
they do not contain sources themselves. By design, this treatment does not include screening 
and evaporation of mini-halos (e.g., Barkana and Loeb (2002)). 

To summarize, the semi-analytic model actually has three phases: the "interior of halos," 
the hot ionized H H regions, and the cold neutral regions. The "interior of halos" phase 
includes all the gas in the universe that is in halos with baryonic mass greater than the 
Jeans mass. We assume that these regions have are virialized with isothermal sphere density 
runs, and therefore include all gas with overdensity greater than Gtt^. In this phase, the 
cooling criterion and the Schmidt law arc used to calculate the ionizing source function, and 
the photoionizations and recombinations arc treated through the escape fraction calculation 
(a Stromgren sphere approximation). Those photons that "escape" this phase then enter 
the "second phase," which are the ionized H H regions. In this phase, a clumping factor 
approach is used, with a high density cutoff at 67r^. The last phase is the cold neutral part of 
the universe, which uses up ionizing photons only when "converted" into the ionized phase. 



(12) 
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2.3. Modeling Gunn-Peterson Absorption 

The statistics of Gunn-Peterson absorption has been reviewed by numerous others 
(Miralda-Escude, Haehnelt, and Rees, 2000; McDonald and Miralda-Escude 2001; Fan 
et al.2002; Songaila and Cowie 2002). To summarize, we begin with the stardard optical 
depth at resonance for a gas with neutral density nni at redshift z: 

^ ^ 3A2p^isAj^nHi(^) .^g. 
87ri/oV^o(l + ^)' + ^^A' 

where A2p--^.is is the decay rate (6.25 x 10^ sec~^), Aq is the Lyman-a wavelength (1.216 x 10~^ 
cm), and the other symbols have their usual meaning in a cosmological context. 

Now, we make the standard assumptions of ionization equilibrium and a uniform ionizing 
intensity in the ionized regions to obtain the neutral fraction. Considering only photoion- 
ization and radiative recombination 

nHir2i J21 = X{X + F/4)/mJ R^Tl^-' pl^\ (14) 

where Y = 0.24 is the helium mass fraction (we assume helium is singly ionized), i?4 a; 
4.2 X 10~^^cm'^ s^^ is the recombination rate at a temperature T = 10"^ K, T^^ = Tq/IO^K is 
the temperature at mean density in units of 10^ Kelvin, and r2i is the photoionization rate 
for J21 = 1. We use the method of Hui and Gnedin (1997) to calculate the temperature at 
mean density after reionization (see also Hui and Haiman 2003). Plugging in values for the 
various other constants (121 = 4.35 x 10"^^ s'^, nip = 1.67 x IQ-^^ g, p;, = 1.88 x 10-^9(1 + 
z)^Qbh'^ g cm~^) gives 

„„,.7.63xlO--(M!)!Zi;!^!(i±i)!em-3 (15) 

or 

^_ 0.316 Q,^fe3(i + ,)6r-0.7 ^^ 
y/Qoil + Z)^ + Qa J2I 

where T4, A, and J21 depend on the redshift z, but only A also depends on spatial location. 
It is helpful to reformulate this for Qo(l + z)^ ^ ^a, and scaled from Q^/i^ = 0.02 and 
Qoh'^ = 0.14: 

r^,.,36(l±±)'\m(9t]lY^^^ (17) 

Vl + 5.5y V Qoh^ \ 0.02 y J21 ^ ^ 

Defining the optical depth for a uniform medium (A = 1) is 

l+z V /01T/'SJJi2\'T4-°' 
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we can define the mean transmitted flux ratio at a given redshift 

/*CXD 

% = (e--"^') = Q / Py(A)e-^"^'ciA, (19) 

JO 

where Q is the volume filling factor for the ionized regions (we have now adopted the standard 
notation — Chin and Ostriker (2000) used "/" as the volume filling factor). Note that strictly 
speaking, there should be a high density cutoff. However, since the exponential has a power 
of A^, for large values of A, the optical depth is very large. The high-density end of the PDF 
contributes very little to the integral, and the results are thus not sensititive to this cutoff. 
In fact, at the redshifts of interest, the integral is usually dominated by values of A < 1, as 
was noted by Barkana (2001) and others. 

Our semi-analytic model provides the values of T4, and J21, and from the previous 
section, we have a prescription for determining Py(A). Thus, we have a complete model 
from which to calculate the expected Gunn-Peterson absorption. 

3. Compciring Model Predictions and Observations 

Having calibrated our model using various hydrodynamic simulations, we are left with 
one free parameter, euv, in addition to the cosmological parameters. Below, we first consider 
the case where euv is constant and derive a constraint on crgilo ^, similar to those from rich 
clusters, and compare the results and predictions to those found by WMAP. We consider 
joint constraints from WMAP and the SDSS quasars on the spectral index n as well as the 
time-depedence of euv- 

3.1. Combining SDSS Quasar Data 

The quasar data we use combine the compilation at ^abs ^ 5.5 Songaila and Cowie 
(2002) (Table 2) and based on measurements of six SDSS quasars at 2; > 5.7 from Becker 
et al. (2001), Fan et al. (2003), and White et al. (2003) (Table 3). Note that we did not 
include the J1044-0125 iz = 5.74) data in Table 3 because it is already incorporated into 
the compilation of Songaila and Cowie (2002). 

In order to combine these data, we must take into account the fact that at redshifts 
z < 5.6, the uncertainties in the transmission data are dominated by intrinsic scatter, while 
at higher redshifts, the uncertainties are dominated by measurement errors. Since our model 
predicts the mean transmission, our likelihood function must properly account for both 
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measurement error as well as the scatter. For a compilation of transmissions, such as that 
by Songaila and Cowie (2002), the contribution of the data Dc — {{zj, Tj)} to the hkelihood 
is simply 

P{%\D,) = exp 

C^|,mean = (^|,scatter/^- (21) 

The contribution to the likelihood from the compiled data is estimated this way. 

For the individual transmission data, especially in the case where both measurement 
error and scatter are important, the simplest way to do this is to estimate scatter CTgcatter and 
add in quadrature with the measurement error a^eas to obtain the total uncertainty in the 
mean for each data point. In particular, for individual transmission data — {{zj, T!,)}with 
known (Gaussian) errors Umeas and (Jscatter, probability of a mean transmission as a function 
of redshift % is 

P{%\Di) = exp 

Note that if the measurement error is neghgible, then Equation 22 reduces to Equation (21). 

To estimate CTgcatter, we use the Songaila and Cowie (2002) compilation for z < 5.6 (interpo- 
lated), and use the data themselves (binned) to determine cTscatter at higher redshifts. These 
results are also tabluated in Tabic 3. The total likelihood function, C, then is given by the 
product of the two separate likelihoods. 



a 



J, mean 



(20) 



j,meas ' j,scatter 



(22) 



3.2. SDSS Queiseir Contraints with a Constant UV-EfRciency Compared with 

WMAP 

For this initial comparison, we fix the other cosmological parameters to their WMAP- 
only best fit values: n = 0.99, Qhh'^ = 0.024, and h = 0.72. We use the WMAP-only results 
from Spergel et al. (2003), Table 1, to ensure clarity in the data underlying our analysis. 
Note, however, that these last two values have additional independent lines of support, as 
noted in Spergel et al. (2003). 

We use a Bayesian method to determine the cosmological constraints. The priors wc 
use are flo e [0.15,0.40] and ag G [0.5, 1.0]. This region essentially bounds the 99% region 
reported for the combined CMB and cluster analysis of Melchiorri et al. (2002). 

In order to constrain Qq and as separately from euv, we marginalize over the latter 
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parameter to yield the marginalized likelihood distribution: 

C{Qo,a8) = J C{ilo,a8,€uy)d€uy. (23) 

The central values and percentile limits are found by integrating over C{flo, ag). 

The result of this analysis is presented in Figure 1, where we plot the likelihood contours 
in the plane flo — as- As can be seem by the figure, there is a considerable degeneracy in 
this plane similar to that derived from rich clusters. Thus, the cosmological constraint from 
Gunn-Peterson absorption can be summarized by 

agl^S'^ = 0.33 ±0.01, (24) 

where the error term is statistical error only (68%). 

Given the complex physics of the reionization model, it is difficult to determine a mean- 
ingful systematic error-bar. However, we did investigate the effects of changing h and n. 
The net result was to shift the right-hand-side of equation (24) so that 

/(] 72\^'^~^^ / n \ -0.85 
(0.33 ± 0.01) (!^) (j^) , (25) 

SO that lower levels of h or n increase the right-hand-side. Thus, lowering h or n slightly 
would make the Gunn-Peterson constrain more consistent. 

These results for a constant UV efficiency are somewhat discordant with WMAP's 
marginalized value of crgQ^-^ = 0.48 ± 0.12, although only at the 1.25-cr level. They are 
consistent with the cluster-determinations of Bahcall et al. (2002) of crgOo^ = 0.33 ± 0.03 
(note different exponent), although again at about the 1-a level. 

The actual value of euv bears some discussion. Although we have left this parameter 
as completely free, there certainly are astrophysical constraints on its value. In our 95% 
region, we find that euv ~ 2 x 10~^ (see left panel of Figure 2). For a Scalo mass function 
with metal-enriched stars (l/20th solar metalicity), the value calculated from population 
synthesis is ~ 5 x 10~^ (Wiyithe and Loeb 2002). Given the uncertainties and the relative 
simplicity of our model, the correspondence is quite remarkable. 

The more striking inconsistency is that in all cases, the optical depth to electron scat- 
tering, 

-1000 



J dz—caTUe, (26) 



where iTt = 6.652 x 10 ^^cm^ is the Thomson cross-section and Ue is the electron density, 
is too low. Our model finds that Tes ~ 0.06 (see right panel of Figure 2), whereas Spergel 
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et al. (2003) and Kogut et al. (2003) report WMAP-only best fit values of Tes = 0.17 ± 0.07 
and = 0.17 ±0.04, respectively. Several articles written in the wake of these results (e.g., 
Haiman and Holder 2003; Ciardi, Fcrrara, and White 2003; Cen 2003) have discussed this 
inconsistency and indicated that star formation must begun much earlier than previously 
thought. We consider this issue in the following analysis. 

3.3. Combined WMAP and SDSS Constraints with a Time- Varying 

UV-EfRciency 

There are good astrophysical reasons to believe that euv may effectively increase with 
redshift. For instance, the first generation of stars would have been metal-free, and stellar 
models predict them to have an significantly higher UV output per baryon (about 4x higher 
was reported by Wiyithe and Loeb (2002)). In addition, the creation of metals and henceforth 
dust, which would obscure ionizing sources, may also lead to lower effective efficiencies at 
lower redshift than at higher redshift. 

Thus we now consider constraints on n and the time-dependence of euv from WMAP 
and SDSS quasars jointly. To limit the dimensionality of the parameter space, we keep the 
baryon abundance and the Hubble constant at their WMAP best fit values, and in addition 
fix Qq^^ = 0.14 to its WMAP best fit value. For each value of n and Tes, WMAP predicts 
a unique best-fit value of the normalization erg (based on the code provided by Verde et al. 
2003 and the accompanying data files from Hinshaw et al. 2003 and Kogut et al. 2003). The 
normalization given approximately by (for our fixed values of h, floh^, and flbh^) 

a^e~^'- = 0.765 + 0.6(n - 1), (27) 

where formula (27) is good to better than 0.5%. Note that at the WMAP best-fit values of 
n — 0.99 and Tes — 0.17, this formula gives the WMAP best-fit value of erg — 0.9. 

We assume the following heuristic form for the time-dependence of the UV-efficiency: 

euv = euv,o(l + ^e-^*/^*-»), (28) 

where /* is the fraction of baryons in stars. Thus, the luminosity of each ionizing source is 
given by 

L{M,) = (1 - /*)M6c2eesce*euv,oC (l + ^e"^*/^--'*) . (29) 

The new parameters A and /*,crit define the time-dependence of the efficiency. The essential 
features are that at very high redshift, the total efficiency is {1 -\- A) times greater than at 
lower redshift, with the transition occuring at /* ~ /*,crit- This form is motivated by the 
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notion that the first stars, being metal-free, should have had a higher effective efficiency, but 
that as metals were produced by this generation of stars and reinjected back into the IGM, 
the efficiency would decrease. However, we emphasize that this model was simply selected 
heuristically, and was not based on any detailed model of star formation feedback on the 
UV-efficiency. For instance, it should be noted that the "enhancement" factor may be a 
combination of changes in the intrinsic efficiency euv and the resolution factor related 
to the fraction of cooling gas that forms stars. To give a sense of the meaning of /*,crit; 
Table 4 presents the transition redshifts -Ztrans where /* ~ (ln2)/* crit so that euv(-2;trans) ~ 
euv,o(l+^/2). 

We fix the value of /*,crit = 10^'^ (we discuss the sensitivity to this parameter below). 
Thus, for each point in the n — Tes plane, our procedure searching parameter space is as 
follows: (1) Pick value of A; (2) Adjust euv,o so that Tes,modei (calculated by model) is 
consistent with Tes; (3) Stop if is minimized, otherwise go back to step (1). For each point 
in the n — r^s plane, we have a value of a best fit X^{n, Tes)- We treat as approximately the 
— 21n£, where jC{n, Tcs\f*,ciit) is the likelihood function, and integrate (with uniform priors) 
to find find the confidence regions (the results are insensitive to the priors). 

The effective optical depth for several models is shown in Figure 3. The three illustrative 
"good-fitting" models are in the 95% region, and the "badly-fitting" model is outside the 
region. Clearly, high values of n and Tes are ruled out by the data. 

Considering all the constraints together, the results are shown in Figure 4. Our analysis 
strongly favors a narrow range (95% errors): 

Tes = 0.11^°:°^ (30) 

n = 0.961°:°! (31) 

for /*,crit = 10~^. This range is consistent with the WMAP results, also shown in the Figure. 
The effect of changing /*,crit; is only to shift the Tes constraint up or down. For instance, the 
constraint for /*,crit = 3 x 10~^ is Tes — ^■^'^-om^ '^^^^ change in the n constraint. The 
results for increasing /*,crit are in the opposite direction, decreasing the constraint on Tes by 
~ 0.02 for a factor of 10 increase in /*,crit- 

4. Discussion 

In order to understand the quasar constraints, let us return to the constant efficiency 
case. Consider the following "subset" of the quasar data: 

Ql Requiring -^overlap ^ 6 and only using the Lya constraint from quasar J1148-I-5251 {zem — 
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6.43); 

Q2 Using the Lj(3 constraint from SDSS 1030+0524 (^em = 6.28); 
Q3 Using the Songaila and Cowie (2002) constraints at z < 5. 

Figure 1 also shows the effects of combining these three constraints. The intersection 
(Q1+Q2+Q3) subset of data appears to adequately depict the upper limit of the combined 
dataset. In particular, these three sets of data contain much of the "information" in the full 
treatment. The lower limit of the combined dataset actually also contained within the 
constraint data in Q3, but changes as a function of -Zoveriap- Since Q3 is the union of all these 
constraints for all values of -^overlap; the lower limit does not appear. 

This illustration shows that given the three degrees for freedom Qq, erg, and euv, the 
quasar data constrains to essentially a line in this space. This is because for each flo, the 
data (Q1+Q2) at 2; ~ 6 and the data (Q3) at 2; ~ 4.5 provide strong, essentially independent, 
constraints on the remaining parameters erg and euy. Of course, the "length" and "width" 
of the the hne (i.e., the exact shape, as illustrated in Figure 1) depend on the hkelihood 
function in detail. 

Now let us return to the time- varying efficiency case. For each value of /*,crit, we have 
four parameters n, t^s, euv,o, and A. Heuristically, we now have three constraints on the 
output: the data (Q1+Q2) at 2; ~ 6, the data (Q3) at 2; ~ 4.5, and the self-consistency of Tes 
(value "in" = value "out"). Therefore, once again, we expect a "line" in the 4-dimensional 
parameter space (for fixed /*,crit)- As before, the "length" and "width" of the the line will 
depend on the likelihood function in detail. 

We note also that the implied values for the normalization in our 95% region is ag ~ 
0.83;^o;q5 (based on WMAP best-fit), which for our assumed value of gives agil^^ — 
0.381q;q25, j^st overlapping with the cluster measurements crgO^-^ — 0.33 ± 0.03. 

The next question is what are the implications as to the implied values of the UV- 
efficiency? Do the efficiencies make sense? Table 5 shows a number of parameter combina- 
tions from the 95% region. The implied efficiency at high redshift is on the order ~ 10"*^; a 
transition occurs at 2; = 15 ~ 20, and the efficiency at low redshift is euv,o = 10~^'^ ^. The 
span of these efliciencies encompasses those typically calculated through population synthe- 
sis of 10~^ Given uncertainties of factors of a few in the gas collapse fraction (our 
resolution factor e*), dust absorption (at lower redshift), etc., the values of the efficiencies do 
not seem unreasonable. They certainly do not approach the upper bound for conversion from 
nuclear reactions of ~ 10~^. Sokasian et al. (2003) and Wiyithe and Loeb (2003) discuss 
further the star formation efficiency as relating to Population II and Population III stars. 



-16- 



Our results are consistent with the results of Cen (2003) in that we find that to reach 
T"es ^ 0.17 requires that the spectral index is positively tilted with n > 1.02. In this case 
the effective UV- efficiency was at least 100 x greater at z 3> 6. However, our calculations 
indicate that models with such high n (and hence high power spectrum normalization erg) 
are inconsistent with quasar transmission measurements at redshift z < 5. 

Like other authors (e.g., Haiman and Holder 2003, Cen 2003), we find that a value of 
= 0.17 is inconsistent with constraints at z < 6 for simple models of reionization, and 
may require more exotic (though not necessarily implausible) methods for creating ionizing 
photons, such as mini-quasars or an X-ray background. However, we conclude that WMAP 
data taken as a whole and quasar observations at 2; < 6 in fact are entirely consistent 
for reasonable values and time-dependence for the UV-efficiency. Our analysis shows the 
importance of taking into account the significant degeneracy between t^s and n determined 
by WMAP. 

Having found a consistent set of evolutionary models, let us describe their properties in 
slightly greater detail. We consider the following two models: a model with constant UV- 
efficiency with h = 0.72, Qo = 0.27, n = 0.99, erg = 0.64, euv = 2.2 x 10-^ and Tes = 0.06); 
and a model with variable efficiency with h — 0.72, = 0.27, n — 0.96, as — 0.827, 
euv = 7.7 X 10-6(1 + 17e-^*/^°"'), and Tes = 0.11. These are both near peak of the likelihood 
distribution for Qq — 0.27. 

The Gunn-Peterson optical depths are already shown in Figure 3. The best-fit models do 
not show much difference, as expected since the parameters were fit to these data. Figure 5 
shows the reionization properties for the two models, although the last phase of reionization 
occurs rapidly in both models, the phase in which the average neutral fraction drops from 
unity to ~ 0.1 takes a significantly longer time in the model with a variable UV-efficiency. 
This is a "necessary" part of the model in order to achieve a higher electron scattering optical 
depth. The transition from the higher to the lower UV-efficiency is clearly seen in the stellar 
baryon fraction /* at redshift z ~ 15. We note that, even in the case of Tes = 0.11, there are 
not two distinct epochs of ionization, in contrast to the calculations of Cen (2002a). Rather, 
there is an extended period during which the ionization increases from 1% to 10% before the 
rapid phase change at 2; ~ 6. 

The thermal properties are show in Figure 6. Again, the rise in the global volume- 
averaged temperature is much more gradual with the variable UV-efficiency. In addition, 
in this case, the final temperature is higher, but this is due to the higher power spectrum 

normalization ((Tg). These temperatures are somewhat lower that the "peak" temperatures 
derived by Hui and Haiman (2003), but this is probably due to the "sudden" reionization 
model used in those calculations. A more gradual reionization transition would smooth out 
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the peak, and would appear consistent with our calculations. For instance, Hui and Haiman 
(2003) consider a "stochastic" reheating process their derived temperatures at z ~ 4 are 
similar to ours. 

Finally, in Figure 7 we show the star formation rate and the electron scattering optical 
depth to redshift z. Star formation rates at ^ ~ 4 from field galaxy measurements have been 
reported at about ~ xq-2.5±o.5 Mpc^^ (Steidel et al. 1999), while a recent lower limit 

of ~ 10-3-^ M0 yr-i Mpc'^ has been reported at 2; ~ 6 (Standway, Bunker, and McMahon 
2003). Our model results imply that there is exists a large population of unobserved sources 
at z ~ 6. We not that quasar observations alone require that Tcs has reached value of ~ 0.05 
at 2; ~ 7. Finally, the electron scattering optical depth shows that in the constant efficiency 
case, the full optical depth is reached by 2; ~ 10, while in the variable efficiency case, the full 
optical depth is not reached until 2; ~ 20. This is an observational signature that could be 
detected by future CMB experiments. 



5. Summary and Conclusions 

We derive constraints on several cosmological parameters based on observations of Gunn- 
Peterson absorption in high-redshift quasars and WMAP observations (Bennett et al. 2003). 
We use a semi-analytic model for reionization (Chiu and Ostriker 2000) that takes into 
account a number of important physical processes both within collapsing halos (e.g., H2 
cooling) and in the intergalactic medium (e.g., H2 cooling, Compton cooling, and photoion- 
ization heating). The model is also calibrated to hydrodynamic simulations. We also develop 
a method for estimating the gas PDF, which is important for properly calculating the mean 
absorption in the IGM, as a function of cosmological parameters. 

We find that the Gunn-Peterson absorption data provide constraints on the power spec- 
trum at small scales in a manner similar to that derived from the cluster mass function. 
Assuming that the efficiency of producing UV photons per baryon is constant, the con- 
straint takes on the form UgfiQ^ = 0.33 ± 0.01 assuming a flat, A-dominatcd universe with 
h = 0.72, n = 0.99, and Vlbh'^ = 0.024. The best fit for the WMAP data (marginalized 
over all parameters) is reported as crgf^o^ = 0.48 ± 0.12, which differs by slightly more than 
1 — a. However, the derived value for the optical depth to last-scattering Tes ~ 0.06 differs 
signficantly from the WMAP-determined value of Tes — 0.17 ± 0.04. 

Since the WMAP constraints on Tes are somewhat degenerate with the value of the 
spectral index n (Spergel et al. 2003), we then let the primordial spectral index n float while 
flxing the best flt WMAP determination of Qq^^ — 0.14 and normahzing the power spectrum 
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using WMAP. In addition, we allow the UV-efficiency to have time-dependence. Combining 
the WMAP constraints with the quasar transmission data, our analysis favors then a model 
with Tes = O.llig:^^, n = OMtom, implying erg = 0.83t[J:{];^ (all at 95% confidence), and an 
effective UV-efficiency that was at least ~ lOx greater at 2; ^ 6. If future observations 
conffrm this range for the optical depth to electron scattering Tcs, then it would appear that 
no more "exotic" sources of UV-photons are necessary. We are unable to find a model that 
is consistent with all observational and physical constraints that has an electron scattering 
optical depth within the 1-a range given by the WMAP team. If the EE data finally require 
a value Tes ^ 0.17, then more exotic sources of early ionizing photons will be required that 
those considered in this paper. 

The authors would like to thank Renyue Cen kindly providing his simulation outputs, 
and Martin Rees, Massimo Ricotti, and David Spergel for useful discussions. XF thanks 
support from the University of Arizona and a Sloan Research Fellowship 

A. Approximate Model for Escape Fraction 

In this appendix, we consider a simple Stromgren sphere approximation for the escape 
fraction. The basic conditions for star forming halos in our semi-analytic model are as 
follows: 

1. In halos greater than the Jeans mass, the baryonic mass will collapse to a singular 
isothermal sphere with = ^A^ph^r^/r)^, where ~ 178 is the virial overdensity, ph 
is the mean baryonic density, and r„ is the virial radius. Here the total baryonic mass 
is Mb = fA^PbT^. 

2. In halos where the cooling time tc is less than the dynamical time td, there will be star 
formation. 

3. The rate of UV photon production per unit baryonic mass is ^ = c^e^euyEQ^f^^ , where 
euv is the mass to UV energy efficiency and Eq — 13.6 eV. Here e* is the fraction of 
cooling baryons in the halo that are forming stars (the "resolution" factor) , and euv is 
the fundamental conversion efficiency from stellar baryons to UV photons. 

In regions that are optically thick, so that all photons are absorbed locally, then local 
ionization equilibrium would imply that Op^ = n^x'^a, where nn is the total hydrogen density, 
X is the ionization fraction, and a is the recombination coefficient. This implies that the 
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ionization fraction is x — \^9pb/a/nii oc r. Therefore, as r ^ 0, the ionized fraction also 
tends to 0. 

To check for self consistency, consider the optical depth over a central region. The 
average density over the sphere r is (p&(r)) = A^p(,(r^/r)^. Optical depth over this region is 
approximately (r) ~ {nn{l — x))ar, where a is an effective cross section. Since (nn) ~ r~^, 
this implies that (r) ~ (1 — x)r~^. Therefore, as long as (1 — x) does not vanish, then the 
optical depth will become very large as r — > 0. Above, we established that a; as r ^ 
for an optically thick region. We therefore have a self consistent picture for the central part 
of the halo. 

Now what is the condition so that the entire halo is optically thick? For a completely 
neutral halo, r ~ Aypi,ary/mcs. Plugging in numbers gives r ~ r„(l+2;)'^/10 kpc (in physical 
units). In comoving units, r oc (1 + z)^. Therefore at any redshift before full reionization 
> 6), all halos of concern are optically thick if fully neutral. They only become optically 
thin if the neutral fraction is small. 

Now consider a model based on a central source approximation. Let S{r) be the number 
of ionizing photons emitted by a central source which pass through a sphere of radius r. The 
standard equation for S{r) is given by 

dr 

Using the definitions above, we obtain 



= -Anr nfiX a. (Al) 



^ = -4. f (A2) 

or \ Srrics ) 

Now consider, the number of photons produced within r, S{r) = 6M{< r). Differentiating 
this with respect to r gives 

dS{r) 4:7^ r. ^ 9 / * ON 

= -eA^p,rl (A3) 

Now let us make the ansatz that all these photons can be considered to be radially emitted. 
Then these two equations can be combined to yield 

Because we know that the halo is optically thick if a; < 1, let us consider the following 

as 

dr 



prescription: as long as a; < 1, all photons are absorbed locally, so = 0. This leads to the 
expression for the ionized fraction as a function of radius 

x{r) — Min 



A^pbOi ' 



(A5) 
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If we define quenching fraction 



then 77 > 1 means that all of the photons are absorbed within the halo. If < 1, then the 
ionized fraction has reached its maximum value of 1 at the "quenching radius" = rjry. 
Note importantly that rj depends only on redshift through the background density pb and 
the dynamical time in 6, and is essentially the same for all halos. (Actually, there is a weak 
dependence on the halo temperature through the recombination rate a.) 

If < 1, then we use equation (A4) with x — 1 and the boundary condition S{rq) — 
to solve for S{r) from to r^. The result is 

S{r,) = —eA^p^rUl - vf = Mb9{l - r^f. (A7) 

Thus, the factor (1 — riY is the fraction of UV photons created in the halo that actually 
escape — the escape fraction Cesc- At high redshift, when > 1, this factor is — all the 
photons are consumed within the halo. 

To illustrate, consider the "best fit" model with time- variable UV-cfficicncy described in 
the main text. Using a temperature of 10^ Kelvin, Vtoh? — 0.14, Vt^h?' — 0.024, the quenching 
fraction is 

77 \ ' ■ (A8) 

In our "best fit" model, we fixed e* = 0.03, while euv ranged from ~ 0.8 x 10"^ at low 
redshift to 1.3 x 10"^ at high redshift. At very high redshift, we have 

This implies that at very high redshift, when > 1, all the photons are absorbed locally so 
that the escape fraction is essentially 0. This is to expected since densities are much higher 
at high redshift. At redshifts near transition ~ 15, when euv ~ 6 x 10"^, we have 

ri ^ [-^j , (AlO) 

implying an escape fraction ~ 0.2, comparable to values used in other semianalytic models. 
At 2; ~ 6, we have 

V « [-^) , (All) 

implying an escape fraction eesc ~ 0.03. These values are all within the range used by others 
in semi-analytic models (e.g., Cen 2003; Wiyithe and Loeb 2002; Haiman and Holder 2003). 
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Fig. 1. — Summary of cosmological constraints in the — ag plane for a constant UV- 
efficiency. The constraints Ql, Q2, and Q3, as described in the text, are labeled in the 
shaded regions ("+" in this context denotes "intersection"). The thick solid lines are the 
formal 68% and 95% contours using all the quasar data. 
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Fig. 2. — The values of the UV-efficiency euv (left panel) and the optical depth to electron 
scattering Tcs (right panel), overlaid on the quasar data constraints (Figure 1), for a constant 
UV-efficiency. 




Fig. 3. — Plot of effective Lyman-o; optical depth Tee = — ln(T) from Songaila and Cowie 
(2002) compilation (boxes, error bars are cr^ean), SDSS quasars (diamonds, errors do not 
include cxscatter), and model predictions for three model runs within the 95% region of the 
n — Tcs plane for /*,crit = 1 x 10~^, and one model run outside of the region (i.e., ruled out 
by the data). The values of (n, Teg), in decreasing likelihood, are thick solid line:(0.96,0.11), 
short dash: (0.93,0.08), long dash: (0.98,0. 13), dot-dash: (0.99,0. 14). For comparison, the thick 
gray line is the best-fit model for a constant efficiency at VIq = 0.27 and n = 0.99 (cxg = 0.64). 
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Fig. 4. — Constraints (68% and 95% contours) from quasar observations for a time-varying 
UV-efficiency in the n — plane, for WMAP-normalized models with /*,crit = 1 x 10~^. 
Also shown (short dashed line) is approximatly the 68% constraint from WMAP (Spergel 
et al.2003, Figure 5). 
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Fig. 5. — Some basic reionization properties of the "best fit" models with a constant UV- 
efficiency (left panel, t^s = 0.06) and with a time-dependent UV-efficiency (right panel, 
Tes = 0.11) as a function of redshift z. Show are the fiUing factor (solid line), neutral hydrogen 
fraction in the ionized region (short dashed) and the global average neutral hydrogen faction 
(long dashed) , and the fraction of baryons in stars /* (long-short dashed) . 
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Fig. 6. — Thermal properties of the "best fit" models with a constant UV-efficiency (left 
panel, t^s — 0.06) and with a time-dependent UV-efficiency (right panel, Tgs = 0.11) as a 
function of redshift z. Shown are the volume- average temperature in the ionized region 
(short dashed) and the global volume-average (long dashed). 
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Fig. 7. — Star formation rate (left panel) and electron scattering optical depth to redshift 
z (right panel) for the constant UV-efficiency (thick gray line, Teg = 0.06) and the variable 
UV-efficiency (thin solid line, Tes = 0.11) "best fit" models. 
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Table 1. Fit Parameters for Gas PDF Derived from Simulations 



Redshift 




Cen'' 






/3 &o 


6 


2.50 


1.09 


2.52 2.6 


5 






2.41 3.0 


4 


2.48 


1.53 


2.25 3.6 


3 


2.35 


1.89 




2 


2.23 


2.54 





''Miralda-Escude et al. 2000 
''Cen 2002b 



Table 2. Compilation of Transmission Data from Songaila and Cowie (2002) 



Redshift 


Mean Transmission T 


(^scatter 


Cmean = Cscatter/x/jV 


N 


4.09 


0.352 


0.352 


0.027 


15 


4.34 


0.334 


0.334 


0.020 


20 


4.61 


0.260 


0.260 


0.017 


15 


4.93 


0.162 


0.162 


0.022 


5 


5.20 


0.107 


0.107 


0.022 


8 


5.ol 


0.071 


0.071 


0,011 


7 
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Table 3. Transmission Data 



Redshift Range 


Transmission T ± (Tmeas 


^'"scatter 


Reference 


SDSS quasar 


Comments 


4.84 - 5.00 


0.1334 ± 0.0011 


0.046 


a 


J0836+0054 




5.00 - 5.17 


0.0809 ± 0.0011 


0.055 


a 


J0836+0054 




5.17-5.33 


0.0523 ± 0.0008 


0.055 


a 


J0836+0054 




5.33 - 5.50 


0.0692 ± 0.0010 


0.029 


a 


J0836+0054 




5.50 - 5.66 


0.1185 ±0.0011 


0.029 


a 


J0836+0054 




5.25 - 5.41 


0.1324 ±0.0036 


0.029 


a 


J1030+0524 




5.41 - 5.58 


0.0996 ± 0.0033 


0.029 


a 


J1030+0524 




5.58 - 5.74 


0.0418 ± 0.0033 


0.020 


a 


J1030+0524 




5.74 - 5.95 


0.0242 ± 0.0038 


0.020 


a 


J1030+0524 




6.0 — 6.17 


r\ r\r\^ r\ I r\ r\r\r\c\ 

0.0010 ± 0.0009 


0.0005 


b 


J 1030+0524 




6.0-6.17 


0.0043 ± 0.0088 




b 


J1030+0524 


Lyman-/3 


5.95 - 6.15 


0.0031 ± 0.0149 


0.0005 


c 


J 1048+4637 




5.50 - 5.70 


0.0541 ± 0.0014 


0.029 


a 


J1148+5251 




5.70 - 5.90 


0.0139 ± 0.0019 


0.020 


a 


J1148+5251 




6.0-6.10 


0.0 ±0.0063 




b 


J1148+5251 




6.0-6.10 


0.0 ± 0.335 




b 


J1148+5251 


Lyman- /3 


6.10 - 6.32 


0.0 ± 0.0021 




b 


J1148+5251 




6.10 - 6.32 


0.0 ±0.051 




b 


J1148+5251 


Lyman-/3 


5.00 - 5.17 


0.1429 ±0.0027 


0.055 


a 


J1306+0356 




5.17-5.33 


0.0922 ± 0.0025 


0.055 


a 


J1306+0356 




5.33 - 5.50 


0.0936 ± 0.0025 


0.029 


a 


J1306+0356 




5.50 - 5.66 


0.0679 ± 0.0027 


0.029 


a 


J1306+0356 




5.66 - 5.83 


0.0699 ± 0.0033 


0.020 


a 


J1306+0356 




5.55 - 5.75 


0.0383 ± 0.0258 


0.020 


c 


J 1630+4027 





Note. — Lyman-a unless otherwise noted 
s-Becker et al. (2001) 

''White et al. (2003). For J1148+5251, these are conservative upper Umits based on the hypothesis of 
Ly-a emission from an intervening galaxy. 

<=Fan et al. (2003) 



Table 4. 2;trans as a Function of n and /*,crit for Semi-analytic Runs 



/* ,crit 














0.94 


0.97 


1.00 


1.03 


1.06 


3 X 10- 


-5 


16.5 


20 


24.5 


30 


35 


1 X 10" 


-4 


14.5 


18 


22 


27 


32 


3 X 10" 


-4 


12.5 


15.5 


19 


23.5 


28 


1 X 10" 


-3 


10.5 


13 


16 


20 


24 



-31 - 



Table 5. Parameter Combinations in the 95% Confidence Region 



"^es /*,crit — 10 ^ 

n euv,o(xlO-*) A euv.high zCxlO"") 



0.09 


0.94 


0.12 


7.8 


0.92 


0.10 


0.95 


0.096 


12 


1.1 


0.11 


0.96 


0.077 


17 


1.3 


0.12 


0.97 


0.058 


26 


1.5 



